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Abstract 

\u 

The  capabilities  of  classical,  minimum  energy,  and  linear  predictive  array  processing 
algorithms  to  estimate  the  bearings  of  two  equal-energy  sources  is  examined.  Signal 

coherence  is  shown  to  affect  adversely  the  resolution  and  detection  capabilities,  as  well 

as  the  bias  characteristics,  of  all  three  algorithms.  For  linear  arrays  of  equally-spaced 
sensors,  the  superior  resolution  capability  of  the  linear  predictive  algorithm  is 
demonstrated.  The  value  of  utilizing  prediction  elements  in  the  center  of  the  array  to 

resolve  very  closely— spaced  source  bearings  is  demonstrated.  However,  the  linear 

predictive  algorithm  is  least  capable  of  detecting  highly  coherent  sources.  A  tradeoff  is 
established  between  resolving  capability  and  sensitivity  to  finite  averaging.  Conditions 
are  established  which  indicate  which  algorithm  is  best  suited  to  anticipated  levels  of 
signal  coherence  and  averaging.  The  estimates  of  source  bearing  produced  by  each 
algorithm  are  shown  to  be  asymptotically  biased.  The  bias  produced  by  the  classical 
beamformer  is  most  severe,  while  the  minimum  energy  beamformer  produces  the  least  bias.  ^ 

i 

!♦  lB*JQ 

The  determination  of  the  bearings  of  distant  sources  of  acoustic  energy  iu  a  noisy 
ocean  environment  is  a  primary  objective  of  a  passive  sonar  system,  A  group  of  acoustic 
sensors  arranged  in  a  known  spatial  pattern  (an  array)  is  deployed  to  record  the  ac£ustic 
field.  Classical  or  Bartlett  beamforming,  minimum  energy  (ME)  adaptive  beamforming  ,  and 
linear  predictive  (LP)  processing  algorithms  are  commonly  employed  to  estimate  source 
bearings  from  array  data  [2],  These  processing  procedures  result  in  beasr-pat terms,  plots 
of  estimated  incident  energy  versus  bear! ng-of-look  (see  figure  1). 

Three  criteria  are  used  to  evaluate  the  performance  of  a  processing  algorithm  with 
respect  to  source  bearing  estimation.  The  first  is  resolution.  When  a  single  propagating 
plane  wave  impinges  on  the  array,  a  global  maximum  occurs  om-tsrget,  i.e.,  at  the  bearing 
of  incidence.  When  two  plane  waves  are  incident,  their  presence  may  or  not  be  evident  in 
the  beam-patterns.  If  separated  in  bearing  by  a  sufficient  amount,  their  bearings  will  be 
resolved,  in  which  case  the  beam-patterns  exhibit  two  distinct  maxima  as  in  figure  1.  On 
the  other  hand,  the  beam-patterns  may  fail  to  resolve  the  source  bearings;  in  this  case, 
the  beasr-patterns  display  a  single  broad  maximum  located  at  some  intermediate  bearing, 
Ihe  beam-patterns  corresponding  to  two  unresolved  sources  can  greatly  resemble  those 

•This  research  was  supported  by  Office  of  Naval  Research  Contract  N00014-81-K-0565. 

*  This  beamforming  algorithm,  due  to  Capon  [1],  is  often  referred  to  ss  the  maximum 
likelihood  method  in  the  literature.  Actually,  the  beam  does  not  maximize  a  likelihood 
function;  thus  the  customary  terminology  is  misleading. 
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corresponding  to  a  single  source.  Resolution  analysis  assesses  the  conditions  nnder  which 
an  array  processing  algorithm  is  ccpable  of  resolving  source  bearings.  A  nsefnl  measure 
of  the  resolution  capability  of  an  array  processing  algorithm  is  the  signal-to-noise  ratio 
it  requires  to  resolve  two  closely-spaced,  equal-energy  sources* 

The  beam-patterns  produced  by  the  processing  algorithms  exhibit  small  local  maxima 
off-target  (figure  1).  These  off-target  ripples  are  referred  to  as  sidelobes.  The  second 
criterion  is  the  degree  to  which  the  sidelobes  are  small,  so  that  they  are  not  confused 
with  peaks  corresponding  to  incident  energy.  Thus,  detection  analysis  assesses  the 
conditions  under  which  the  number  of  targets  present  in  the  acoustic  field  can  be 
determined  accurately. 

Finally,  a  processing  algorithm  can  both  detect  and  resolve  sources  but  yield 
inaccurate  estimates  of  the  bearing  of  each  target.  Bias  analysis  assesses  how  much  the 
estimates  of  source  bearing  deviate  from  the  true  target  bearings. 

Cox  [3]  analyzed  the  capability  of  the  Bartlett  and  HE  array  processing  algorithms  to 
resolve  two  equal-energy,  temporally  nncorrelated,  sources.  He  demonstrated  that  the  ME 
algorithm  can  resolve  closely-spaced  source  bearings  at  a  lower  signal-to-noise  ratio  than 
can  the  Bartlett  algorithm.  Studying  the  single  source  case,  Seligson  [4]  observed  that 
the  performance  of  the  ME  algorithm  is  quite  sensitive  to  deviations  of  the  signal 
wavefront  from  its  assumed  planar  shape.  For  example,  the  superposition  of  temporally 
correlated  plane  waves,  which  occurs  in  multipath  environments,  produces  a  net  wavefront 
which  is  not  planar.  As  multipath  propagation  is  prevalent  in  the  ocean  [S],  there  is  a 
need  to  understand  how  source  correlation  affects  the  resolving  capabilities  of  array 
processing  algorithms.  This  paper  pursues  an  analysis  similar  to  that  of  Cox  to  study  the 
capabilities  of  Bartlett,  ME,  and  LP  algorithms  to  detect  and  resolve  coherent  (i.e. 
correlated,  narrowband)  signals  using  linear  arrays  of  equally- spaced  sensors  (LES 
arrays).  Furthermore,  bias  of  the  bearing  estimates  produced  by  the  three  algorithms  is 
studied. 

Our  detect  ion/ re  solution  and  bias  analyses  are  based  on  the  assumption  that  an 
infinite  amount  of  ergodic  array  data  is  available,  thereby  allowing  the  cross  spectral 
correlation  matrix  to  be  determined  exactly  by  averaging  over  an  infinite  period  of  time. 
Consequently,  these  results  are  not  sufficient  to  assess  the  performance  capability  of  an 
algorithm  used  in  practical  situations.  To  provide  some  insight  into  the  utility  of  each 
processing  algorithm,  guidelines  are  established  which  indicate  the  amount  of  averaging 
necessary  to  suppress  variability  of  the  beam-patterns  due  to  finite  averaging. 

2.  Signal  Model 

The  acoustic  field  present  in  an  ocean  environment  can  be  modelled  as  a  finite 
superposition  of  P  signal  plane  waves  propagating  with  different  bearings  Jt  and  energies 

added  to  a  background  noise  field  having  zero  mean.  Assuming  that  the  m**1  sensor  is 
located  at  position  its  output,  sampled  with  period  T,  is 

x  <nT)  -  n  <nT)  +  >  s  (nT  -  «j  )  m  -  0,1,2 . M-l  (la) 

*  “  p£i  P  c-p  m 

M  is  the  number  of  sensors  in  the  array,  c  is  the  speed  of  propagation,  and  J[  is  a 
dimensionless  unit  vector  pointing  along  the  bearing  of  propagation.  We  restrict  our 
attention  to  array  geometries  which  exhibit  symmetry  through  the  origin,  i.e. 
ju_1  ■  -|  ,  as  in  figure  2.  By  definition,  the  signals  seen  at  the  origin  of  the  array 
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have  zero  delay.  Such  a  model  is  said  to  be  phase  centered.  Evaluating  the  discrete 
Fonrier  transform  of  (la)  yields 


J(f) 


o0(f)N(f)  + 


9  (f)S  (f). 
P  P 


(lb) 


where  each  component  in  the  complex  M-dimensional  vectors  corresponds  to  an  array  element, 
<rA(f)N(f)  is  the  vector  of  discrete  Fonrier  transforms  of  the  noise,  S  (f)  is  the 

0  -  p 

f 

direction  vector  corresponding  to  sonrce  p  with  components  S  *  exp(-j2n  k  • z  ),  and 

pn  c  p  # 

2 

9  ( f)  is  the  discrete  Fonrier  transform  of  s  (nT).  on(f)  is  the  variance  of  the  noise  at 
p  p  o 

I 

frequency  f  and  N(f)  is  normalized  so  that  E(N  N)  =  M, 

The  following  simple  narrowband  model  is  used  to  describe  the  signals: 


j(2nf  nT  + 

s  (nT)  *  or  e  P 
P  P 


where  o  is  a  real  deterministic  amplitude  constant,  f  is  the  deterministic  center 
P  P 

frequency,  and  p  is  a  random  phase  variable.  The  signal  DFTs,  evaluated  at  the  analysis 
P 

frequency  f,  are 

j[pp  +  n(fp-f)T(N-l)] 

9  (f)  -  o  e 
P  P 

where 

_  sin(n(fn-f)TN) 

°p  "  °p  sin(n( f° -f)T)  *  (2> 

P 


Referring  to  figure  2,  the  components  of  the  source  direction  vectors  can  be  written  as 

S  «  exp(ja  ),  Array  geometry  determines  the  parameters 
pm  p,m 


a 


P,m 


1  'J 
P 


nr 

uf  «»•<»,  -  v 


wh,r,  r  to,  y  in  tb,  polar  coordinate,  of  the  a”  sanaor  aod  9  ia  tba  bearing  of  tbo 

■  B  p 

pih  source.  X  is  the  common  wavelength  of  the  narrowband  signals. 

2 

Eaoh  processing  algorithm  involves  the  evaluation  of  a  functional  : 

2 

For  notational  convenience,  the  dependence  of  the  beam  energies  on  frequency  is 
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(3a) 


i  ,  ,  rx 

Pm  *  l£„E  I*  •  (3c> 

LPq  I  q  I 

K  is  the  cross  spectral  correlation  natriz  [2].  J  is  the  directios-of-look  rector  witn 

elements  I  -  ezp(ja  )  where  a  »  -j2n  ^  ,  and  k  is  the  bearing-of-look.  U  is  the 

m  m  ■  c  ■  Q, 

prediction  element  rector  containing  a  one  in  the  q**1  element  and  zeros  elsewhere.  The 
usual  value  of  q  for  a  LES  array  is  0,  i.e.  the  prediction  element  is  at  the  end  of  the 
array;  here  we  consider  the  variation  of  using  other  prediction  elements.  typically.  a 
power  of  -2  rather  than  -1  is  used  for  the  LP  algorithm.  Here,  nse  of  the  -1  power 
ensures  that  the  units  of  (3a.b.c)  agree  and  allows  direct  comparison  of  the  resolving 
capabilities  of  the  algorithms  [2].  An  equivalent  formula  for  the  LP  beam  energy  which  is 
more  convenient  for  our  analyses  is 


PLP  (i> 
q 


(3d) 


-1  -1  '  -1 
where  -  B  x. 

The  foundation  upon  which  this  theoretical  study  rests  is  the  assumed  form  of  the 
cross  spectral  correlation  matrix. 

*  -  e[h’] 


Using  equation  (lb),  one  obtains 


K 


The  cross  correlations  of  the  signal  DFTs  are 


suppressed.  Ve  examine  the  resolution  and  detection  capabilities  at  s  single  analysi 
frequency. 
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L  p  ®  j  p  ■  L  J 


-f  )T(N-1) 


The  signal  coherence  coefficients  are  defined  to  he  the  normalized  cross  correlations 

Ek<) 


c  -  - 

pm  a  a 
pm 


and  are  analogous  to  correlation  coefficients.  It  is  easily  shown  that  the  magnitude  of 

* 

each  coherence  coefficient  lies  between  zero  and  one.  Furthermore,  C  =  C  ,  When 

pm  mp 

signals  p  and  a  differ  by  only  a  constant  deterministic  phase  and  amplitude,  they  are 

perfectly  coherent,  i.e.  |C  |  *  1,  This  situation  can  arise,  for  example,  when  multiple 

pm 

propagation  paths  (with  fixed  differences  in  path  length)  exist  between  a  single  source 
and  the  array. 

We  consider  the  case  where  the  noise  is  spatially  white,  has  zero  mean,  and  is 
independent  of  the  signals.  Subject  to  these  assumptions,  the  cross  spectral  correlation 
matrix  can  be  written  as 

B  «  a20I  +  SCS'  (5) 

th. 

where  the  PxP  coherence  matrix  C  is  composed  of  elements  C  ;  S  is  an  MxP  matrix  whose  p 

pm 

column  is  ,  Using  a  well-known  matrix  inverse  formula  [6],  the  correlation  matrix 
inverse  can  be  expressed  as 


! 

9 

-1 

R-1  -  \  I  -  S 

2  2 

9 

CS  . 

°0  °0 

•o  J 

Here,  we  are  interested  in  the  two  source  case.  For  convenience,  the  signal 
coherence  C12  is  expressed  as  *  c  “  Icle^.  The  array  aigaal-to-moiae  ratios  (ASNK)  of 

2  2  2  2 

the  two  sources  are  defined  to  be  ■  M(o^/oq)  and  ■  M (o^/oq),  respectively.  The 

geometric  mean  of  the  ASNRs  is  “  j Aj •  *n  this  case,  the  correlation  matrix  and  its 
inverse  can  be  expressed,  respectively,  as 


.\W .  > 


L'i:-*. rU  iit  i  on/ 

A.-aliability  Co<j 
;Avall  aDd/o^ 
Dist  I  Special  j 
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*■  °0  1  +  itw'i  ♦  *  «**a]  +  r*A 


m  M  r 

1  ^1  »  T.2  -  '  * 

2  1  “  M -1^1  “  M  ICJS152  +  *  -2~1  " 
°0 


*2  . 


The  similarity  in  form  of  the  correlation  matrix  and  its  inverse  is  readily  apparent  from 
equations  (7)  and  (8).  The  parameters  A^,  and  c  are  given  by 

|l  +  A,,(l-Ic|2)  I  A, 

\  *  - = - - - r  (9a) 

l+A1+A2+2A12Re  0*008^  ,^2)  +AJ  A2  <1- J  c  I2 )  1-cos2  (j^.Sj) 


+  Ajd-lcl2^ 

l+A1+A2+2 A12Re  |c*cos  < ^ »J2) | +AJ Aj (1- I c |2 )  |l-cos2  (j^,  Sj) 

c  -  A^d-lcl^cosCj^,  Sj) 

0  *  1 Z  •  1  1  — » 

xj(l  +  A1(l-Ic|2))(l  +  A2  < 1- 1 c 1 2 ) ) 


A^  and  A2  are  the  modified  ASNRa  of  the  two  sources.  The  geometric  mean  of  the  modified 

ASNRs  ia  A12  •  \|AjA2>  As  with  coherence,  the  magnitude  of  the  modified  coheremce,  c.  lies 
between  zero  and  one.  The  modified  coherence  can  be  expressed  in  terms  of  its  magnitude 

**  |**|  j0  w  a# 

and  phase:  e  **  |c|e  .  For  equal-energy  sources  A  »  ■>  Aj  -  Aj2  and  A  *  A^  ■  A2  »  A^. 

Finally,  for  the  two  aource  case,  the  matrix  B  2  appearing  in  equation  (3d)  can  be 
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written  as 
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*  [  Jai  a  »  _J®t  a  *1  ^2  ^°2  a  *  ~^®2  a  *1 

v»  «l*  ’ Vi ♦  •  'VJ-ttI*  V* +  •  vj 

^12  ^°2,a  '  ~  “ja2,a  '  ~  ^°1  q  '  ~*  , q„  *1 

-  « lc  •  Vi  *  ••  V, ♦  ••  V2  ■* c  •  VJ 

m  m  r  m 

M  MM  M  M 

^1  *  ^12  <1  >  -le  '  ^2  ' 

+  j^Ml  +  f[*  *1*2  +  •  ^  -2-1 J  +  ^2-S2 


MM  M 

where  the  parameters  A^,  and  p  are  given  by 


~  j  -ja,  -ja.  !2 

A  .  Il  •  1,4  +  A  7.  2,ql 
A1  |A1®  A12Ce  | 


(11a) 


~  ~  _ja2.fl  -  “j0l,q| 


(lib) 


\\~  "J«i-  ~  ~-jo2,qIf~  Ja2.q  ~  ~*J°l,q| 

Le|lAle  +A12Ce  JlV  +A12ce  j 


(11c) 


M  M 

and  are  the  donblj^-modif  led  ASNSa  of  the  two  sonrceso  The  geometric  mean  of  the 

M  I  M  M  MM  M  M 

SV  I  tfW  iw  0^0  0^0 

donbly-nodif ied  ASNRs  is  A^  ■  \|AlA2*  ^or  e,lual"energy  sources  A  ■  A^  *  A^  *  A^. 

1*  Pyttvtlw  JUUl  RptPlPtlPP  Si£  §SSLS£  Bearing 

To  study  the  effect  of  source  bearing,  signal-to-noise  ratio,  and  source  coherence  on 
the  detection  and  resolution  capabilities  of  the  various  array  processing  algorithms,  we 
restrict  our  attention  to  the  case  where  two  equal-energy  signals  impinge  on  the  array. 
For  convenience,  we  define  the  complex  cosine  between  two  complex  vectors  to  be 

ii 

=  Hill  Hill’ 
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1.1.  Resorption 

When  the  sources  are  resolved,  the  beam  energy  evaluated  at  either  target  bearing 
must  be  larger  than  the  beam  energy  evaluated  between  the  target  bearings.  The  be  tween- 
targets  direction-of-look  vector,  W^,  is  defined  to  have  elements 

Wgn  =  expjj(a^  m+a2  m)/2j,  *?or  *  ^ES  array,  when  the  source  bearings  are  symmetric  about 

array  broadside,  the  between-targets  bearing-of-look  lies  halfway  between  the  target 
bearings,  i.e.  at  broadside.  The  criterion  used  for  resolution  of  the  sources  is  that  the 

2 

ratio  of  on-target  to  between-target  beam  energies  exceed  a  threshold  value  of  it  /8.  This 
quantity,  known  as  the  Rayleigh  resolution  limit,  is  the  ratio  obtained  by  the  Bartlett 
algorithm  using  a  LES  array  composed  of  many  elements  when  the  signal-to-noi se  ratio  is 
large  and  the  source  bearing  separation  is  equal  to  one  beam-width.  A  beam-width  is  the 
bearing  separation  between  the  center  of  the  main  lobe  and  the  first  null  of  a  classical 
beaai-pattern. 

The  following  formulae  exactly  express  the  dependence  of  the  on- to-be tween- target 
beam  energy  ratios  on  ASNR,  coherence,  and  bearing  separation. 


_  ,  .  i  +  aI 

PBART(°n)  1 

111  +  cos2 (5  , S  )  + 

L  1  1 

2  Ic  lcos()»)cos(S1,S2) 

(12) 

PBART(betWeen>  ' 

1  +  2Acos2(Wq, S^) 

[l  +  lclcos(p)] 

pm!od) 


PjjE^be  tween) 


1  -  2Acos2(J0#J^) [l  +  lclcos(p)] 

-  a|i  +  cos^Sj.Jj)  +  2  Ic  IcosijOcos^,  S^)  j 


(13) 


PLP  (°n) 

_ q. .  , 

PLP  (between) 

q 


-  4A 


1  +  2A<1  +  cos(w))cos  (J0.51) 
- ;]  + 


cos 


a„  “®,  1 

l .  o 


2  J 


1  +  a|i  +  cos^Sj.jSj)  +  2cos(?)cos(S1,S2)  j 
2A  lc|(cos(a  -a  -e)+cos(«))  +  1  +  cos(a„  -a,  )cos(J.  ,,$«)  I 

L  2*<i  2»q  1  2  J 


1/2 


(14) 


These  expressions  are  valid  for  all  phase-centered  arrays  which  exhibit  symmetry  through 
the  origin.  F'  *  such  g*  aetries,  the  complex  cosines  between  the  source  and  direction- 
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of- look  vectors  are  real  valued.  The  effects  of  array  geometry,  number  of  sensors,  and 
source  bearing  separation  are  contained  entirely  in  these  cosine  terms. 

The  complexity  of  expressions  (12,  13,  and  14)  prevents  a  clear  understanding  of  the 
effects  of  ASNR,  coherence,  and  bearing  separation  on  resolving  capability.  To  gain  such 
an  understanding,  numerical  solutions  were  obtained  for  the  resolvent  array  signal-to- 
noise  ratio,  the  minimum  ASNR  which  results  in  resolved  sources,  A  LES  array  was  assumed, 
and  the  sources  were  assumed  to  be  located  symmetrically  about  array  broadside.  The 
resolvent  ASNR  is  a  function  of  the  magnitude  and  phase  of  the  coherence  as  well  as  source 
bearing  separation. 

Figure  3a  compares  the  minimum  ASNR  necessary  for  a  ten  element  array  to  resolve  the 
bearings  of  incoherent  sources  using  the  Bartlett,  HE,  and  LPq  processing  algorithms.  Of 

the  three,  the  LPq  algorithm  requires  the  least  ASNR  to  resolve  source  bearings.  The  ME 

and  LPq  algorithms  are  both  capable  of  resolving  arbitrarily  closely- spaced  source 

bearings  if  the  ASNR  is  high  enough.  In  contrast,  the  Bartlett  algorithm  is  incapable  of 
resolving  bearings  spaced  more  closely  than  one  beam-width  regardless  of  the  ASNR. 
Figures  la,  b,  and  c  illustrate  Bartlett,  ME,  and  LPq  beam-patterns,  respectively,  which 

just  resolve  the  indicated  incoherent  sources.  In  each  case  the  ASNR  is  lOdB  and  the 
bearing  separation  is  the  smallest  resolvable  (from  figure  3a).  A  sensible  way  to  compare 
the  capabilities  of  the  algorithms  to  resolve  coherent  sources  is  to  select  the  coherence 
phase  which  maximizes  the  resolvent  ASNR  for  each  source  bearing  separation.  Thus  the 
worst-case  resolvent  ASNR  is  sufficient  to  guarantee  bearing  resolution  regardless  of 
coherence  phase;  sources  with  particular  coherence  phases  can  usually  be  resolved  at  a 
somewhat  lower  ASNR.  Figures  3b  and  c  compare  the  worst-case  resolvent  ASNR  for  the  three 
algorithms  when  the  coherence  magnitude  is  ,7  and  .99,  respectively.  As  the  magnitude  of 
the  coherence  increases,  the  resolvent  ASNR  for  the  ME  and  LPQ  algorithms  increases.  The 

LPq  algorithm  retains  its  advantage  over  ME  and  Bartlett  processing  with  regard  to  bearing 

resolution.  The  LPq  algorithm  consistently  requires  roughly  15dB  less  ASNR  than  the  ME 

algorithm  to  resolve  closely-spaced  source  bearings.  The  Bartlett  algorithm  again 
exhibits  a  minimum  separation  below  which  source  bearings  cannot  necessarily  be  resolved; 
this  limit  increases  as  the  magnitude  of  the  coherence  increases. 

As  discussed  by  Johnson  [2],  the  linear  predictive  algorithm  offers  flexibility  in 
the  choice  of  prediction  element  q.  Figure  4  compares  the  minimum  ASNR  necessary  for  a 
ten  element  array  to  resolve  the  bearings  of  incoherent  sources  using  the  linear 

predictive  algorithm  with  prediction  elements  0  to  4.  As  a  function  of  prediction  element 
q,  the  resolvent  ASNR  is  symmetric  about  the  center  of  the  array.  For  closely-spaced 

source  bearings  (6^4°)  the  center  prediction  elements  (q=4,5)  require  the  least  ASNR  to 
resolve  the  source  bearings.  Figures  5a  and  b  show  the  LP^  and  LPq  beam-patterns, 

respectively,  which  result  when  incoherent  sources  separated  by  2°  impinge  on  a  ten 

element  array;  the  ASNR  is  21  dB.  As  expected  (figure  4),  the  LP.  algorithm  resolves  the 

4 

source  bearings  whereas  the  LPq  algorithm  does  not.  More  widely-separated  bearings 

(4°<6£22°)  are  most  easily  resolved  by  prediction  elements  at  the  ends  of  the  array 

(q*0,9).  The  sharp  peaks  in  the  plots  for  prediction  elements  3-6  correspond  to 

situations  in  which  these  LP  algorithms  produce  spurious  peaks  in  their  beam-patterns  at 
bearings  between  the  actual  target  bearings.  Figure  6  illustrates  this  phenomenon  when 
q*3.  While  the  source  bearings  are  visibly  resolved,  peak  bias  and  the  spurious 
intermediate  peak  combine  to  confound  the  resolution  criterion.  For  all  prediction 
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elements,  as  source  coherence  increases  in  magnitude  from  zero  to  one,  the  resolvent  ASNR 
rises  quite  uniformly  as  a  function  of  hearing  separation.  For  a  ten  element  array,  the 
resolvent  ASNR  for  sources  having  .7  coherence  magnitude  is  6dB  greater  than  for 
incoherent  sou  ces;  the  resolvent  ASNR  for  sources  having  .99  coherence  magnitude  is  15dB 
greater. 


In  addition  to  resolving  source  bearings,  it  is  necessary  for  a  processing  algorithm 
to  detect  the  presence  of  sources  by  exhibiting  peaks  in  the  beam-pattern  which  stand  out 
against  off-target  ripple.  The  criterion  used  for  detection  of  the  sources  is  that  the 
ratio  of  on-target  to  highest-sidelobe  beam  energy  exceed  a  threshold  value  of  two.  This 
ratio  is  established  in  two  stages.  First,  the  ratio  of  on-target  beam  energy  to  noise 
beam  energy  is  established.  By  definition,  the  noise  beam  energy  is  evaluated  when  no 
sources  are  present,  i.e.  A^  =  A^  =  0.  Second,  a  relationship  is  established  between  the 

highest-sidelobe  and  noise  beam  energies.  Together,  these  relationships  are  used  to 
assess  each  algorithm's  detection  capability. 

The  following  formulae  exactly  express  the  dependence  of  the  on- targe t- to-noi se  beam 
energy  ratios  on  ASNR,  coherence,  and  bearing  separation. 


PBtliT(°°) 

PBART<“i“ 


y  ■  1  +  *[,  +  +  2  Ic  lco»(,)cos(S^,  S^) 


Pm(on) 


P^Cnoise) 


-a- 

1  -  A  1  +  cos2^,^)  +  2lclcos(?)cos(S1,S2) 


PLP  (0n) 
_ 9 _ 

PLP  *noise) 
R 


1  +  All  +  cos  +  2cos(?)cos(S1,  S2> 


i 

I  -  2a[ I c  I  (cos (< 


L  L 


■«  -a,  -p)+cos(d))  +  1  +  cos 

m#  ^  1 0  Q 


)cos(S  ,S  ) 
i,  q  l,  q  l  zj 


These  expressions  are  valid  for  all  phase-centered  arrays  which  exhibit  symmetry  through 
the  origin.  The  effects  of  array  geometry,  number  of  sensors,  and  source  bearing 
separation  are  contained  entirely  in  the  cosine  terms. 

For  LES  arrays  it  is  straightforward  to  show  that  the  worst-case  (largest)  ratios  of 
highest-sidelobe-to-noise  beam  energy  are: 


-10- 


De  6mf  and  Johnson 


Bearing  Estimation 


PBART^hl*he,t 

16A 

(18) 

PB*ST(“oi*') 

*  1  + 

** 

PME(highest  sidelobe) 

1 

P^noise) 

!  _  ±$2 

9n2 

(19) 

,  (highest  sidelobe) 
a 

l 

PLP  *noise) 
q 

l 

JL 
“  3jt 

mt 

A  +  -A 
3n. 

(20) 

Equations  (15,  16,  and  17)  in  conjunction  with  (18,  19,  and  20)  establish  lower  bounds  on 
the  on-target-to-highest-sidelobe  beam  energy  ratios  for  LES  arrays.  Osing  these  bounds, 
numerical  solutions  were  obtained  for  the  detectable  array  signal-to-noi se  ratio,  the 
minimum  ASNR  which  causes  the  ratio  of  on-target  to  highest- si  delobe  beam  energies  to 
exceed  the  threshold  value  of  two.  As  with  the  resolvent  ASNR,  the  detectable  ASNR  is  a 
function  of  the  coherence  magnitude  and  phase  as  well  as  source  bearing  separation.  The 
detectable  ASNR  of  the  Bartlett  algorithm  is  virtually  unaffected  by  increasing  source 
coherence.  The  LP  algorithm  requires  the  least  ASNR  to  detect  incoherent  sources,  while 
the  Bartlett  algorithm  requires  the  highest.  For  coherent  sources,  the  worst-case 
coherence  phase  maximizes  the  detectable  ASNR.  Thus  the  worst-case  detectable  ASNR  is 
sufficient  to  guarantee  source  detection  regardless  of  coherence  phase;  sources  with 
particular  coherence  phases  can  usually  be  detected  at  a  somewhat  lower  ASNR.  As  source 
coherence  increases,  the  worst-case  detectable  ASNR  for  the  ME  and  LP  algorithms 
increases.  For  highly  coherent  sources,  the  LP  algorithm  requires  the  highest  ASNR  to 
guarantee  source  detection,  and  the  Bartlett  algorithm  the  least.  All  of  the  algorithms 
require  a  higher  SNR  to  detect  widely-separated  sources  than  to  resolve  them.  Conversely, 
all  of  the  algorithms  require  a  higher  SNR  to  resolve  closely-spaced  sources  than  to 
detect  them.  Figures  7a,  b,  and  c  show  the  minimum  ASNR  required  by  the  Bartlett,  ME,  and 
LPq  processing  algorithms  to  guarantee  both  detection  and  resolution  of  equal-energy 

sources  with  varying  degrees  of  coherence. 

4-  Bill 

Even  when  targets  can  be  detected  and  resolved  in  a  background  of  white  noise,  the 
corresponding  peaks  in  the  beam-patterns  may  not  occur  at  precisely  the  true  source 
bearings  (i.e.,  they  may  be  biased).  As  passive  sonar  systems  are  often  used  to  localize 
targets  at  great  ranges,  small  errors  in  bearing  estimates  can  lead  to  large  errors  in 
position  estimates.  In  addition,  non-zero  bias  would  cause  the  resolvent  and  detectable 
ASNR  of  section  3  to  be  somewhat  higher  than  necessary  for  visible  peak  resolution  and/or 
detection.  When  a  single  source  is  present,  all  three  algorithms  yield  an  unbiased 
estimate  of  the  source  bearing  when  the  correlation  matrix  is  known;  when  multiple  sources 
are  present,  the  estimates  are  generally  biased.  In  this  section,  formulae  are  presented 
which  approximately  describe  the  bias  induced  in  one  target  peak  by  a  second  target  for 
the  Bartlett,  1C,  and  LP  beam-patterns.  We  emphasize  that  the  bias  results  are  based  on 
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the  assumption  that  the  cross  spectral  correlation  matrix  is  known  exactly.  Thus  the 
formulae  indicate  the  asymptotic  (with  respect  to  averaging)  bias  inherent  in  each 
processing  algorithm. 

We  restrict  onr  study  of  bias  to  LES  arrays.  The  components  of  the  direction-of-look 
vector  are  V  *  exp(-j (M-l-2m)a/2) ,  and  the  components  of  the  source  direction  vectors  are 


lm 


exp(-j(M-l-2m)aj/2)  and  S2b 


exp (-j(M-l-2m) 0^/2) .  The  transform  baarimg-of— look,  a. 


is  related  to  the  angular  bearing-of-look,  6,  by 

u  =  n  sin€, 

where  d  is  the  sensor  separation  and  X  is  the  signal  wavelength.  The  source  transform  and 
angular  bearings  are  similarly  related. 

The  transform  bias  of  the  peak  corresponding  to  source  1  is  defined  to  be 

A  A 

Acij  =  a?  -  a^,  where  is  the  estimate  of  the  transform  bearing  of  source  1,  i.e.  the 

location  of  the  actual  beam-pattern  peak.  The  following  notation  will  be  used  throughout 
the  discussion: 

sin(M(a-a, ) /2) 


8l(a)  *  "  tin( (a-a^)/2) 


sin(H(a-a2>/2) 

,2<«>  -  Hco.[i.32]  -  ,u«.-.2>/2) 


Using  equations  (7)  and  (8),  the  Bartlett  and  ME  beam-patterns  can  be  expressed  as 
functions  of  transform  bearing  as  follows: 


PBART(ct) 


jr  rn 


J0  [M  +  ^gl(a)  +  2^lc|cos(^)gl(o)g2(a)  +  Y*2(a>] 


(21) 


i-l 


[PME(o>]  *  *  I  1  "I  “  “T 


M 


-  "j^gjta)  “  2-~|c  lcos(5»)g1(o)g2(a)  -  ~g2(a) 


(22) 


When  the  bias  is  small  compared  to  a  beam-width,  2n/M,  Taylor  series  expansions  of  g^  g2> 
2  2 

*1'  *2'  *0<*  *i®2  t0  Sflcoa^  order  about  will  be  accurate  in  the  neighborhood  of  the 

maximum.  Substitution  of  these  expansions  into  (21)  and  (22)  and  setting  the  derivatives 
with  respect  to  u  equal  to  zero  at  the  maximum  yields 
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Bartlett:  Aa, 


[A12lclcoa(#)M  +  A2*2(oi,]«2(ai) 
A1Mg1(o1>  +  +  *2(al)g2<0l>] 

+  A12lclcoa(#)|Mg2(o1)  +  «2(al)gl(al)] 


ME :  Ao^  «* 


^A12lc|cot(0)M  +  A2g2(a2)  g^ctj) 

XlMgl(al>  +  +  g2(0i)g2<0l)] 

+  A^lcIcoslp^Mgj^)  +  82(al)gl(al)] 


Taylor  series  analysis  analogous  to  that  already  described  is  applied  to  the 
quadratic  fora  (3d)  to  obtain  the  LP  bias  expression.  For  notational  convenience,  define 
Cj  (a)  ■  cos(Q(a-a2 ) )  and  S2<a)  ■  sin(Q(a-a0 ) ) .  where  Q  *  (M-l-2q)/2. 


Do  6ruf  ud  Johiioi 


Bearing  Estimation 


LP  :  Aa. 

q  1 


w  m  I 

A^costJOM+Ajg^a^jg'u^ 

[C2<al)82(ol)"QS2(ol)g2(ol)] 
<?)[,'  (•1)-»s2(.1,] 
-  MA12lc|sin(#)[*OC2(a1)-Qg2(o1)] 


-  MA 


-  MA^IcIco* 


Xl*Vai>  +  X2[82(al>82(al)+822(al)] 

+  A12coa(?)rMg2(o1)+gi<a1)g2(a1)j 

~  r «  2 

-  MA1lg1(a1)-MQZ 

"  tf2[C2(ol)8i,(al)-2QS2(ol)82(ol)-Q2c2(°l)82(al)] 
MA12lc|co.(?)|c2(o1)g2(o1)-MQ2C2(o1)+g2(«1>-<l2g2t«1)] 

-  MA12lc|«in(?)|^S2(o1)gJ(a1)-IIQ2S2(o1)-2Qg2(o1)j 


(25) 


Equations  (23,  24,  and  25)  approxiaate  the  effects  of  source  bearing  separation, 
coherence,  signal-to-noi se  ratio,  relative  source  strength,  and  number  of  array  elements 
on  peak  bias.  While  the  expressions  are  complicated,  several  important  observations  can 
be  made  pertaining  to  each  algorithm.  For  the  Bartlett  algorithm,  bias  does  not  depend  on 
SNR.  When  the  sources  are  incoherent,  the  bias  of  peak  1  increases  as  the  strength  of 
source  2  increases.  The  direction,  positive  or  negative,  of  the  bias  depends  on  the 
source  separation.  While  there  are  some  bearing  separations  for  which  the  bias  is  zero, 
generally  it  is  non-zero.  The  behavior  of  ME  peak  bias  as  a  function  of  relative  source 

strength  and  coherence  is  clearly  more  complicated  than  is  the  Bartlett  bias  behavior 

(refer  to  eqs.  (9)  and  (ID).  However,  the  similarity  of  the  ME  and  Bartlett  bias 

expressions  is  obvious.  For  equal  energy  sources  A^  -  Aj,  and  the  only  difference  between 

equations  (23)  and  (24)  is  the  substitution  of  c  for  c.  The  modified  coherence  (c)  of 
sources  which  are  not  perfectly  coherent  tends  asymptotically  toward 

-cos^.jgj)  *  as  A  SNR  increases.  Consequently  the  ME  bias  for  equal-energy 

sources  tends  to  zero  as  ASNR  increases.  However,  for  finite  ASNR  the  ME  bias  is 
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generally  non-zero.  Similarly,  it  can  be  shown  that  the  LP  bias  for  equal-energy  sources 
tends  to  zero  as  ASNR  increases;  for  finite  AS NR  the  LP  bias  is  generally  non-zero* 

The  complexity  of  the  bias  expressions  necessitates  computer  evaluation.  Figure  8a 
compares  the  Bartlett,  ME,  and  LPQ  bias  of  two  incoherent,  equal-energy  sources  using  a 

ten  element  array;  the  ASNR  is  lOdB.  The  approximately  linear  region  of  each  plot  for 
small  source  separations  reflects  the  inability  of  the  algorithms  to  resolve  the  sources* 
At  larger  source  separations,  the  bias  exhibits  oscillatory  behavior.  The  MB  algorithm 
has  the  smallest  amplitude  of  oscillation,  followed  by  the  LPq  and  Bartlett  algorithms. 

In  figure  8b,  the  ASNR  has  been  raised  to  30  dB.  The  bias  of  both  the  ME  and  LP^ 

algorithms  has  decreased  for  all  source  separations;  more  closely-spaced  sources  are 
resolved  than  before.  For  widely-spaced  sources,  the  ME  bias  remains  smaller  than  the  LP^ 

bias.  The  Bartlett  bias  is  unchanged  by  the  increase  in  ASNR.  Figures  9a,  b,  and  c 
compare  the  Bartlett,  ME,  and  LPq  bias  when  the  coherence  magnitude  is  .7  and  the  ASNR  is 

10  dB.  Each  plot  corresponds  to  a  different  coherence  phase.  Clearly  coherence  phase  as 
well  as  source  separation  determine  the  direction  and  magnitude  of  the  peak  bias.  For 
widely-separated  sources,  the  Bartlett  algorithm  continues  to  exhibit  the  largest  bias, 
and  the  ME  algorithm  the  smallest. 

Figure  10  compares  the  LP  bias  of  incoherent  source  bearing  produced  by  prediction 

<1 

elements  0-4  using  a  ten  element  array;  the  ASNR  is  lOdB.  Note  that  the  plot  for  the 
center  prediction  element  ( q=4 )  departs  from  linearity  at  a  smaller  source  separation  than 
the  others,  reflecting  its  superior  resolution  capability.  However,  the  center  prediction 
element  can  produce  larger  bias  than  the  others  for  more  widely-spaced  sources. 

When  the  sources  are  separated  by  at  least  an  amount  corresponding  to  the  first  zeros 
of  (23,  'd  and  25),  the  biases  predicted  by  these  expressions  are  very  accurate.  However 
expressiuu^  (23,  24,  and  25)  do  not  predict  the  biases  of  more  closely-spaced  sources 
accurately.  The  bias  for  unresolved  sources  can  be  a  significant  fraction  of  one  beam- 
width.  Consequently  the  validity  of  the  bias  expressions  for  closely- spaced  sources  is 
suspect  because  they  are  based  on  Taylor  expansions  of  the  beam-patterns  to  second  order. 
Expansion  of  the  beam-patterns  to  third  order  produces  significantly  more  complicated  bias 
expressions  which  do  not  provide  a  correspondingly  great  improvement  in  accuracy. 

1*  MfS9tf  Si  FifllU  Ayyygg^ng 

Results  presented  in  sections  3  and  4  were  based  on  use  of  the  true  correlation 
matrix.  The  effect  of  empirical  computation  of  the  correlation  matrix  by  averaging  the 
available  sensor  output  data  must  also  be  considered.  Assume  that  T=NK  samples  of  the 
sensor  outputs  are  available.  The  data  is  broken  into  K  consecutive  non-overlapping 

segments  each  of  length  N.  The  ktl  data  segment  is  transformed  yielding  a  vector  of 

Fourier  transforms  ,J^(f).  Following  common  practice  [7,8],  the  correlation  matrix 
estimate  is  given  by: 


£ 


k«0 


The  statistics  of  this  random  matrix  have  been  studied  by  Goodman[7],  It  is  important  to 
note  that  £  will  not  be  positive  definite  or  iinvertible  if  E  <  M  because  it  will  have  M-K 
eigenvalues  equal  to  zero. 
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Con aid ar  the  following  narrowband  signal  model: 

_  j(2nf  nT  +  p  (nT)) 

.  (nT)  -*  e  p  p 

P  P 


where  f  is  the  center  frequency  of  source  p  and  0  (t)  is  a  slowly  varying  random  phase. 
P  P 

It  is  assumed  that  p  (t)  varies  slowly  enough  that  it  can  be  considered  essentially 

P 

constant  over  a  length  N  DFT. 


4*1. 


SaksssMa 

Evaluating  the  DFT  of  the  k**1  signal  segment  at  the  analysis  frequency  f  yields 


(f) 


j[*P.k 


+  k2nf  IN  +  r(f  -f)T(N-l)l 
P  P  J 


o  e 
P 


where  #  .  «  e  (kNT)  and  o  is  as  defined  in  equation  (2).  The  coherence  coefficient 

p»  k  rP  P 

between  sources  p  and  m  is  defined  to  be 

C  .  •  <2«> 

pm  ~  0  a  * 

pm 


As  and  are  random  variables,  so  are  the  coherence  coefficients.  Comparing 

p  m 

equations  (26)  and  (4),  it  is  apparent  that  the  expectation  operator  has  been  replaced  by 
a  time  average  of  the  sample  functions. 

One  important  mechanism  which  results  in  source  coherence  is  multipsth  propagation. 
To  see  how  the  coherence  due  to  multipath  depends  on  the  time-bandwidth  product,  K, 
consider  the  two  source  case.  The  effect  of  multipath  propagation  is  modelled  by  writing 


*l.k  "  #2.k 


*  * 


where  p  is  a  constant  random  phase  shift  due  to  the  average  difference  in  path  lengths, 
and  is  a  zero  mean  random  phase  varisble  which  accounts  for  a  time  varying  difference 

in  path  lengths  due  to  platform  motion,  source  motion,  or  motion  of  the  scattering  surface 
(see  figure  11).  For  mathematical  convenience,  assume  that  p^  and  p^  are  statistically 

independent  when  i  p  k,  and  that  the  probability  density  of  p ^  is  uniform  over  the 

interval  l-pn,pn].  The  conditional  expected  value  of  the  coherence  given  p  is 


E(C12I*> 


pn 


sin(n(f1-f2)lNK) 
Ksin(jr(f  -f,)IN) 


j[*  +  nif^mNK-l)] 


(27) 


De  Graaf  aad  Johnson 


Bearing  Estinatioa 


The  nagnitnde  of  the  conditional  expected  value  of  C12  i  a  always  less  than  or  equal  to 

one.  It  tends  toward  zero,  but  not  monotonically,  as  the  time-bandwidth  product  and/or 
the  inter-signal  frequency  difference  and/or  the  randomness  of  the  inter-signal  phase 
increases  (p — >1).  The  coherence  between  sources  which  have  nearly  the  same  center 
frequencies  or  phases  is  reduced  by  increasing  the  amount  of  averaging  K.  Jf  the  signals 
share  the  same  center  frequencies  and  have  a  constant  phase  relationship  then  the 
coherence  magnitude  will  be  one  regardless  of  the  amount  of  averaging  performed.  Even 
when  the  inter-signal  phase  sequence  is  completely  random  (p  =  1)  the  variance  of  the 
magnitude  of  the  coherence  is  1/K,  not  zero.  Thus  sources  which  are  incoherent  in  the 
expected  value  sense  of  equation  (4)  may  not  be  incoherent  in  the  time  average  sense  of 
equation  (26). 

£.2,.  Randomness  of  Beam- Patterns 

Capon  and  Goodman[8]  have  studied  the  statistics  of  the  Bartlett  and  ME  beam-patterns 
when  finite  averaging  is  performed  to  obtain  the  correlation  matrix  estimate,  and 
Baggeroer[9]  has  studied  the  statistics  of  the  LP  spectral  estimator  for  time  series  data. 
These  analyses  are  rigorous  mathematically,  but  fail  to  provide  a  simple  guide  to  the 
amount  of  averaging  required  to  reduce  beam-pattern  randomness  to  acceptable  levels.  By 
analyzing  the  variability  of  on-target  and  noise  beam  energies  with  respect  to  the  non¬ 
zero  signal-noise  cross  terms  in  the  single  source  case  [10],  the  following  guidelines  can 
be  established  which  dictate  the  amount  of  averaging  necessary  for  each  algorithm  to 
produce  beam-patterns  with  a  small  random  component.  These  results  are  valid  for  LEM 
arrays. 


Bartlett : 


(28) 


ME: 


LP: 


K 

I 


>> 

>> 


M  and  K  >>  M2 ~ 
a0 


(29) 

(30) 


Figures  12a,  b,  and  c  show  several  Bartlett  beam-patterns  for  a  ten  element  array 
when  K*=10,  100,  and  1000,  respectively;  the  SNR  is  one.  The  variability  of  the  Bartlett 
beam-patterns  decreases  as  K  increases.  In  all  three  cases  the  amount  of  averaging  is 
such  that  the  criterion  K>>.1  is  satisfied.  For  K=10  the  degree  of  variability  is  greater 
than  might  be  expected;  with  such  small  time-bandwidth  products  the  assumptions  leading  to 
condition  (28)  are  not  entirely  valid;  nevertheless  the  variability  for  K=10  is  not  so 
severe  as  to  obscure  the  target.  The  variability  of  both  the  K=100  and  K*1000  beam- 
patterna  is  small.  Figures  13a,  b,  and  c  show  the  ME  beam-patterns  for  the  same  data. 
The  variability  of  the  ME  beam-patterns  decreases  as  K  increases.  The  K=10  case  does  not 
satisfy  condition  (29),  and  the  variability  is  so  great  as  to  obscure  the  broadside 
target. 

Figures  14a,  b,  and  c  show  several  LPq  beam-patterns  for  a  ten  element  array  when 

K*>10,  100,  and  1000,  respectively;  the  SNR  is  one.  The  variability  of  the  patterns 
decreases  as  the  time-bandwidth  product  increases.  It  has  been  shown  [10]  that  finite 
averaging  can  also  produce  spurious  noise  induced  peaks  off-target  in  LP  beam-patterns; 
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inch  peak*  are  apparent  for  the  K=10  case.  Only  the  K*1000  beam-pattern  has  small 

variability,  as  expected,  since  M^IOO.  The  second  condition  of  (3o)  indicates  that  as  the 
SNR  increases,  more  averaging  mnst  be  performed  to  maintain  a  constant  level  of  noise 
induced  variability  in  the  beam  energies.  Fignre  15  shows  several  LP^  beam-patterns  for  a 

ten  element  array  when  Ks1000;  the  array  SNR  is  20dB.  Comparison  of  figures  14c  and  15 
confirms  that  the  variability  of  the  on-target  beam  energy  increases  as  the  SNR  increases. 
For  large  arrays,  the  second  condition  (30)  is  dominant  because  of  its  proportionality  to 

M2. 

£.  Conclusions 

The  bearing  estimation  capabilities  of  the  beamforming  and  linear  predictive  array 
processing  algorithms  discussed  here  are  adversely  affected  by  source  coherence  and  by 
limited  averaging  in  the  computation  of  the  correlation  matrix.  Furthermore,  the 
performance  of  each  algorithm  is  affected  to  a  different  degree  by  these  conditions.  For 
LES  arrays,  no  one  algorithm  is  superior  to  the  others  in  all  situations.  Consequently, 
the  coherence  and  amount  of  averaging  possible  in  particular  applications  greatly 
influence  the  choice  of  an  "optimal"  array  processing  algorithm. 

The  incremental  effect  of  a  change  in  coherence  magnitude  on  the  detection  and 
resolution  capabilities  and  bias  of  the  ME  and  LP  algorithms  is  greatest  when  the 
coherence  magnitude  is  large  (near  one);  for  the  Bartlett  algorithm,  the  incremental 
effect  is  greatest  when  the  coherence  magnitude  is  small.  Thus,  the  adaptive  algorithms 
(ME  and  LP)  are  less  sensitive  to  low  levels  of  signal  coherence  than  is  conventional 
beamforming.  As  the  magnitude  of  the  coherence  increases  from  zero  to  one,  the  resolving 
capability  of  all  of  the  algorithms  diminishes.  Except  for  perfectly  coherent  sources, 
the  LP  processing  algorithm  is  uniformly  the  most  capable  of  resolving  closely-spaced 
sources,  followed  by  ME  and  Bartlett  processing.  The  superior  resolution  capability  of 
the  LP  algorithm,  compared  with  the  ME  algorithm,  is  somewhat  offset  by  the  greater  bias 
of  its  bearing  estimates.  In  terms  of  statistical  averaging,  the  Bartlett,  HE,  and  LP 
algorithms  do  not  yield  asymptotically  unbiased  estimates  of  source  bearing.  Therefore, 
Craaer-Rao  lower  bounds  on  the  variance  of  these  bearing  estimates  are  suspect  as  they 
usually  assume  that  the  estimates  are  asymptotically  unbiased.  The  version  of  this  bound 
that  pertains  to  biased  estimators  contains  an  additional  factor  of  one  plus  the 
derivative  of  the  bias  with  respect  to  the  parameter  being  estimated  [11].  Even  for  small 
biases,  the  magnitude  of  this  derivative  can  be  large  with  respect  to  one  because  of  the 
oscillatory  nature  of  the  bias.  For  LES  arrays  the  maximum  magnitude  of  this  slope 
appears  to  be  independent  of  the  number  of  array  elements  when  the  ASNR  remains  constant. 
The  ME  algorithm  produces  the  least  bias  in  its  bearing  estimates,  and  that  bias  decreases 
as  SNR  increases.  LP  bias  is  greater  than  ME  bias,  but  also  diminishes  with  increasing 
SNR.  Bartlett  bias  is  the  most  severe  and  is  independent  of  SNR.  The  bias  produced  by 
all  three  algorithms  becomes  more  severe  as  source  coherence  increases.  The  LP  algorithm 
has  a  slightly  greater  capability  to  detect  incoherent  sources  than  the  ME  and  Bartlett 
algorithms.  Coherence  has  negligible  effect  on  the  detection  capability  of  Bartlett 
processing.  In  contrast,  as  source  coherence  increases,  the  detection  capability  of  ME 
processing  decreases,  as  Seligson  [5]  observed.  As  source  coherence  increases,  the  LP 
algorithm  becomes  the  least  capable  of  detecting  sources,  and  the  Bartlett  algorithm  most 
capable.  Thus  a  tradeoff  exists  between  the  resolution  and  detection  of  highly  coherent 
sources. 

The  variation  to  the  LP  algorithm  of  utilizing  different  prediction  elements  proves 
to  be  of  some  value.  Moving  the  prediction  element  from  the  end  to  the  center  of  a  LES 
array  enhances  the  capability  of  the  linear  predictive  algorithm  to  resolve  extremely 
closely-spaced  source  bearings.  However,  the  beam-patterns  produced  by  the  center 
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prediction  elenents  are  nore  likely  to  exhibit  spurious  peaks  when  the  sources  are  less 
closely  spaced.  In  addition,  bias  of  the  bearing  estimates  can  be  greatest  when  the 
center  prediction  elements  are  utilized.  Choice  of  prediction  element  has  no  effect  on 
tne  capability  of  the  LP  algorithm  to  detect  widely-separated  sources.  The  effect  of 
increasing  source  coherence  is  the  same  for  all  prediction  elements. 

In  section  $  constraints  on  the  time-bandwidth  product  were  presented  which  ensure 
minimal  sensitivity  of  the  beast-pat  terns  to  imperfections  in  the  correlation  matrix  that 
result  from  finite  averaging.  The  Bartlett  beam-pattern  is  least  sensitive  to  finite 
averaging  while  LP  beam-patterns  are  most  sensitive.  The  sensitivity  of  the  ME  beam- 
pattern  to  finite  averaging  lies  between  these  extremes.  Thus,  a  tradeoff  exists  between 
resolving  capability  and  sensitivity  to  finite  averaging. 

The  LP  algorithms  are  most  capable  of  resolving  closely-spaced  sources;  however, 
their  high  sensitivity  to  finite  averaging  restricts  their  application  to  environments 
where  large  amounts  of  averaging  are  possible.  The  Bartlett  algorithm  is  least  sensitive 
to  finite  averaging  and  requires  the  least  SNR  to  detect  widely-separated  sources; 
however,  its  resolving  capability  is  very  poor.  The  ME  algorithm  lies  between  the  LP  and 
Bartlett  algorithms  in  terms  of  resolution  and  detection  capabilities  as  well  as 
sensitivity  to  finite  averaging.  In  addition,  the  ME  algorithm  produces  the  least  biased 
estimates  of  source  bearing.  Thus  the  ME  algorithm  seems  best  suited  to  applications 
where  the  amount  of  averaging  possible  is  small  and  the  capability  to  resolve  closely- 
spaced  source  bearings  is  requisite.  The  LP  algorithms  are  of  little  value  in  these 
situations  because  of  their  extreme  sensitivity  to  noise.  The  Bartlett  algorithm  is  also 
of  little  value  here  because  of  its  extremely  poor  resolving  capability.  Table  1  crudely 
suggests  which  array  processing  algorithm  is  best  suited  to  each  of  four  combinations  of 
source  separation  and  amount  of  averaging  possible  when  a  LES  array  is  used. 


Time-bandwidth 

1  Bearing  separation  of  sources  1 

product 

1  Closely  spaced  . 

1 

Widely  separated  1 

1 

Small 

1  | 

®  ME  for  all  j 

1 

Bartlett  for  1 

•  coherences  s 

!  i 

all  coherences  i 

1 

i 

LP  for  all  | 

1 

Bartlett  for 

all  coherences  J 

i  i 

j  coherences 

1 

1  1 

LP  or  ME  for 

Large 

I  1 

I  LP  for 

small  coherences  { 

center  | 

1 

I 

i  i 

j  all  coherences  | 

1 

when  signals  are  j 

1 

L. 

|  very  close  } 

_ 1 

Table  1. 

Suggested  processing  algorithms  for  LES  arrays. 
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ILilStt  Captions 

Figure  1 :  Bartlett  (a),  ME  (b),  and  LPq  (c)  beam-patterns  generated  using  a  ten 

element  LES  array  with  half-wavelength  sensor  spacing.  The  ASNR  is  lOdB,  and  the 
indicated  sources  are  just  resolved. 

Fiaure  2:  Array  geometry  which  is  symmetric  through  the  origin. 

Figure  :  Vorst-case  resolvent  ASNR  for  the  ten  element  LES  array  when  the  coherence 
magnitude  is  0  (a),  .7  (b),  and  .99  (c). 

Figure  £:  Dependence  of  linear  predictive  resolvent  ASNR  on  prediction  element  using 
the  ten  element  LES  array*  The  sources  are  incoherent. 

Figure  £:  LP^  (a)  and  LPq  (b)  beam-patterns  that  result  when  incoherent  plane  waves 

separated  by  2°  impinge  on  the  ten  element  LES  array.  The  ASNR  is  21  dB. 

Figure  £:  LP^  beam-pattern  that  results  when  incoherent  plane  waves  separated  by 

14.023°  impinge  on  the  ten  element  LES  array.  The  ASNR  is  30  dB.  The  peak  at  0°  is 
spurious. 

Figure  Vorst-case  resolvent  and  detectable  ASNR  for  the  ten  element  LES  array  when 
the  coherence  magnitude  is  0  (a).  .7  (b),  and  .99  (c). 

Figure  8:  Bias  of  incoherent  source  bearing  estimates  using  the  ten  element  LES  array 
when  the  ASNR  is  lOdB  (a)  and  30dB  (b) 

Figure  9:  Bartlett  (a),  IE  (b),  and  LP^  (c)  bearing  estimate  biases  for  sources  witn 

coherence  magnitude  .7  using  the  ten  element  LES  array.  Each  plot  corresponds  to  a 
different  coherence  phase.  The  ASNR  is  10  dB. 

Figure  Dependence  of  linear  predictive  bearing  estimate  bias  on  prediction 

element  using  the  ten  element  LES  array.  The  sources  are  incoherent,  and  the  ASNR  i  s  10 
dB. 

Figure  ljL :  Multipath  propagation  in  the  ocean  due  to  specular  reflection.  The 
difference  in  path  lengths  varies  as  a  function  of  time  because  of  the  dynamic  ocean 
surface. 

Figure  l^:  Variability  of  Bartlett  beam-patterns  induced  by  finite  averaging.  Each 
plot  shows  four  beam-patterns  corresponding  to  distinct  estimates  of  the  correlation 
matrix  using  the  ten  element  LES  array.  The  time-bandwidth  product  is  10  (a),  100  (b). 
and  1000  (c).  The  ASNR  is  10  dB. 

Figure  j.3 :  Variability  of  ME  beaai-patterns  induced  by  finite  averaging.  Each  plot 
shows  four  beam-patterns  corresponding  to  distinct  estimates  of  the  correlation  matrix 

using  the  ten  element  LES  array.  The  time-bandwidth  product  is  10  (a),  100  (b),  and  1000 

(c).  The  ASNR  is  10  dB. 

Figure  14 :  Variability  of  LP^  beam-patterns  induced  by  finite  averaging.  Each  plot 

shows  four  beam-patterns  corresponding  to  distinct  estimates  of  the  correlation  matrix 

using  the  ten  element  LES  array.  The  time-bandwidth  product  is  10  (a),  100  (b),  and  1000 

(c).  The  ASNR  is  10  dB. 

Figure  Variability  of  LPq  beam-patterns  when  the  time-bandwidth  product  is  1000 

and  the  ASNR  is  20  dB.  The  variability  of  the  on-target  beam  energy  increases  as  the  SNR 
increa  sea. 
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